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INTRC DUCT ION 


Progress  has  been  made  this  year  on  two  distinct 
aspects  of  the  research  project;  the  use  of  two-dimensional 
spectral  analysis,  and  the  calculation  of  vertical  derivatives 
of  a  surface,  starting  in  each  case  from  an  altitude  matrix. 

The  results  are  presented  in  modular  form,  with  several  reports  on 
each  aspect  of  the  research.  The  first  report  deals  with 
calibration  of  a  spectral  program,  using  artificial  waveforms 
of  known  properties.  Experience  thus  gained  was  used  in 
re-rv:nning  2-D  spectra  for  altitude  matrices  in  the  Digital 
Terrain  Library,  the  subject  of  the  second  report.  Thirdly, 
a  summary  of  the  value  of  2-D  spectra  in  general  geomorphometry 
is  offered.  Two  reports  on  vertical  derivatives  follow;  the 
first  presents  preliminary  results  for  the  Digital  Terrain 
Library,  while  the  second  summarises  derivatives  for  different 
(photo- interpretation)  terrain  classes  within  a  matrix  area. 
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CAblDFATION  TEST 3  ox 


FOR  2-0  S:';X7.  nXk:  YaTS 


Initial  results  from  using  the  2-D  spectral  analysis  program  of  Raynor 
C  1*V7Z) 

and  McCaldcn^ suggested  a  need  for  greater  experience  m  selecting  control 
parameters,  to  which  the  results  arc  quite  sensitive.  The  possibility  of 
bias  in  the  2-D  spectral  domain  (Evan  s  and  i>ain,  197 > i  cd 

desirability  of  obtaining  results  for  known  inputs,  such  as  cosine  wave 
fluctuations.  These  trials  firstly  provide  a  test  of  program  accuracy,  ana 
secondly  suggest  some  approximate  guidelines  for  control  parameters.  The 
outcome  is  to  give  the  Rayner-McCalden  program  a  remarkably  clean  bill  of 
health. 

The  control  parameters  involved  are  G,  the  number  of  points  smoothed 

('tapered')  at  the  edges  of  the  data;  (X-D) ,  the  number  of  zeros  added  to 

2 

each  row  or  column;  and  M  ,  the  number  of  spectral  estimated  required.  D 
is  the  width  of  the  data  set  (or  the  depth;  we  will  deal  here  cr.ly  with 
square  matrices),  which  is  generally  given,  while  N  is  the  width  of  t'.-.s 
working  array  including  the  zero  fill. 

The  degrees  of  freedom  of  an  elementary  spectral  estimate  are  proportio 
to  the  square  of  (D-G )/x  (Rayner,  1971,  197(a).  Increasing  G  or  X  decreases 
the  degrees  of  freedom  per  elementary  estimate,  which  is  undesirable  because 
it  decreases  the  stability  of  such  estimates.  On  the  other  hand,  lack  of 
edge  smoothing  and  zero  fill  produces  serious  'side  lobes’,  spurious  peaks  i 
the  spectral  domain  at  harmonics  of  major  peaks.  Some  compromise  must 
therefore  be  made.  Tukey  (1961)  suggested  a  w'ide  range  of  possible  values 
for  G. 

Similarly,  increasing  the  number  of  elementary  estimates  grouped  into  a 
spectral  band  reduces  the  variance  of  spectral  estimates,  but  increases  band¬ 
width,  i.e.  decreases  resolution  (Jenkins,  1961).  This  'smudging'  may  bo 
reduced  only  at  the  expense  of  spectral  stability.  In  t..e  present  ex -;e rimer, 
we  test  the  effect  of  varying  G  and  (X-D),  but  not  M,  which  is  held  constant 
at  lp.  The  effects  of  varying  two  properties  of  the  generated  surface,  the 
wavelength  and  orientation  of  fluctuations,  are  also  assessed.  In  each  c^se 
the  other  variables  are  held  constant. 

Hence  experiments  1  and  2  are  concerned  with  choice  ,,f  parameters  for 
running  2-D  spectral  analyses,  while  3,  4  and  5  measure  the  accuracy  of  the 
urogram  in  locating  osci 1  aliens  of  known  properties. 

All  these  experiments  utilised  'dummy'  in~.ut ,  fo>-  which  the  results  cou 
be  predicted.  A  short  program  was  written  to  create  a  100  x  100  data  matrix 


containing  a  regular  cosine  fluctuation  parallel  to  one  axis,  r.r-~v'.s:i  :*r. 
made  for  alteration  of  wavelength  and  amplitude.  by  addition  of  a  d.as.. 
constant  the  program  could  rotate  the  cosine  wave.  It  was  also  pass!. .la  to 
create  more  complex  surfaces  by  the  addition  of  two  or  ...ore  wave  co:.  ;.-.tnt;5, 

Wave length  is  expressed  in  data  spacing  units,  i.e.  a  IGG-ur.it  wave  goes 
through  one  cycle  across  the  matrix. 

1 .  Tapering  fcdqe  smoothing) 

\vith  constant  window  width  (  (N-D)  =  30  ),  variance  (<s50)  ,  wave  length  and 
orientation  of  input,  the  extent  of  tapering  was  varied  over  the  interval 
suggested  by  Tukey,  i.e.  from  V,'  to  25"'  of  the  data  set  length.  Tukey's 
cosine  bell  function  was  used  throughout. 

Fig.  1  shows  the  percentage  of  variance  remaining  the  frequency  domain, 
plotted  against  G  expressed  as  a  of  0.  The  graph  shows  an  almost 
linear  decline.  As  it  stands,  Fig.  1  represents  the  steeper  portion  of  the 
graph  of  variance  remaining  against  G.  Kith  a  homogeneous  data  set  the 
variance  lost  depends  directly  on  the  number  of  points  affected  by  ta  .-ring . 

G  rows  and  columns  of  the  data  set  are  tapered  on  each  side,  so  each  increase 
in  G  affects  shorter  rows  and  columns.  Therefore  each  increase  ir.  C  affects 
fewer  additional  points,  and  the  rate  of  loss  of  variance  should  also 
decline:  but  this  is  not  apparent  for  G  less  than  20"'. 

Somewhat  different  relationships  would  occur  if  variation  in  tne  data  sat 
were  distributed  unevenly.  For  example,  if  the  edgso  of  the  dat-  set 
exhibited  relatively  more  variability  than  the  centra!  portions  then  the- 
decline,  would  be  more  rapid  over  lower  values  of  C.  Such  a  situation 
violates  second-order  stationarity,  and  would  bo  difficult  to  remove  by 
preliminary  de-trending . 

A  measure  of  efficiency  of  the  choice  of  parameters  is  given  by  the  amount 
of  variance  allocated  to  the  correct  class  in  the  frequency  domain.  Fig.  2 
shows  percentage  variance  in  the  correct  class  plotted  as  a  function  of  G. 

A  maximum  of  93’55'>''  occurs  when  G=  17'  of  D,  when  just  ur.-or  Y'~‘-  of  the 
original  variance  remains  (with  no  tapering,  50"'.  remains).  However,  it  is 
unlikely  that  the  best  value  of  (G/D)  will  be  constant  for  ,  complex 

data  sets.  It  is  therefore  necessary  to  try  several  different  values,  even 
though  in  a  less  controlled  situation  it  may  be  very  difficult  to  decide 
which  result  is  'best1. 
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Addition  of  Zeros 

With  constant  tapering  (G=10),  variance  v^pO),  wavelength  ar.d  orientation  of 


innut,  the  number  of  zeros  added  to  increase  the  data  set  length  was  altered. 
Not  all  values  of  N  can  be  assessed,  since  the  program  requires  that 
(N/2M)  should  be  an  odd  integer  greater  than  1,  where  M  is  the  number  of 
estimates  required.  It  was  found  that  with  N'=100,  i.e.  with  no  a  ros 
added,  the  variance  in  the  correct  class  was  83. By';  with  N=i-  ')  the 
variance  was  Variance  remaining  in  the  frequency  domain  was 

62.56'\'>  and  31  •  58'’'  respectively,  but  the  'loss  of  variance'  due  to  addition 
of  further  zeros  around  the  margins  of  a  matrix  is  not  as  detrimental  as 
the  loss  due  to  tapering.  No  information  is  lost:  total  variability  remains 
the  same,  but  since  this  is  distributed  over  a  larger  number  of  points, 
variance  is  reduced.  Hence  the  larger  value  of  N  is  the  more  successful. 

3.  Wavelength 

With  constant  window  width  (  (N-D)  =  30  ),  tapering  (G=10) ,  variance  (CrpO) 
and  orientation  the  wavelength  of  the  input  was  varied  over  a  range  that 
extended  from  2.16  to  200  units.  This  demonstrated  the  program's  ability 
to  locate  different  frequencies  precisely  in  the  spectral  domain,  in  both 
1  and  2  dimensions.  It  was  particularly  useful  to  see  the  effects  of 
fluctuations  of  wavelength  greater  than  or.e  quarter  of  the  data  set  length, 
a  condition  analogous  to  that  frequently  met  when  the  stationarity  assumption 
is  violated. 

Initially  a  set  of  decreasing  wavelengths  was  fed  into  the  program;  the 
results  are  given  below. 


Approximate 

wavelength 

(units) 

Total 

variance  in 

data  set 

%  Total 

variar.ee  repairing 
after  edge  smoothing 

r'  variance  ir.  2 

v^o ...a «.  r* 

correct  cell 

17.00 

49.80 

46.37 

92.34 

8.68 

50.02 

37.57 

93-49 

7.00 

49.57 

46.44 

92.53 

5.20 

49.51 

57.73 

93.52 

4 . 34 

49.84 

46. 16 

92.52 

3.70 

49.96 

37.39 

93.50 

3.25 

49.50 

46. 51 

92.57 

2.90 

45.62 

37.67 

93-50 

2.60 

50.07 

46.00 

92.55 

2.36 

50.19 

37.20 

9  3.45 

2.16 

48.99 

47.81 

92.12 

92.12 


Differences  in  total  variance  are  dee  10  truncation  of  the  r.  .  — 

different  heights.  That  is,  the  wave  always  1  bee  ins'  at  a  maxi  mu  t  hut 
its  final  value  depends  upon  the  relationship  of  the  wavelength  to  the 
length  of  the  data  set. 

ICach  wavelength  chosen  coincides  with  a  class  midpoint  and  it  is  clear 
from  column  four  that  the  program  works  well  in  placing  c)2-<jh''  of  the 
variance  in  the  spectrum  into  the  correct  class.  At  least  a  further 
5">  of  the  variance  is  concentrated  in  the  eight  cells  immediately  adjacent 
to  the  correct  one.  This  sometimes  leaves  less  than  2"'  of  the  variance 
distributed  over  the  remaining  355  cells  in  the  frequency  domain  specified 
here,  so  that  'leakage'  from  the  correct  wavelength  and  orientation  dc  .s 
not  appear  serious.  The  shaded  maps  of  the-  spectral  domain,  however,  convey 
a  more  pessimistic  impression.  This  is  because  the:  r  geometric  cias.: 
interval  emphasises  contrasts  which  involve  very  small  proportions  of 
variance  (less  than  O.Olf'  per  cell).  I  a  these  simple  examples  the  Rayr.er- 
McCalden  program  exhibits  considerable  accuracy.  In  each  cast,  the  1-0 
averaged  vector  spectra  peaks  are  quite  sharp,  and  are  located  accurately. 

A  similar  experiment  was  performed  at  higher  resolution  to  see  what  happens 
to  fluctuations  whose  wavelength  falls  near  a  class  boundary. 

The  wavelength  was  altered  in  increments  of  0.36  units  over  a  range  extending 
between  two  class  mid-points  in  the  frequency  domain  (6.3  and  3.07  units). 
The  results  are  given  in  rig.  3  where  the  two  linos  or.  tie  graph  represent 
the  differing  percentages  of  variance  in  each  class  as  the  wavelength  a.ecrs. 
It  is  notable  that  the  decline  away  from  each  class  mid-point  is  not  linear. 
The  values  remain  fairly  constant  until  close  to  the  class  boundary ,  vn.sro 
rapid  change  occurs.  Serious  leakage  of  variance  into  adjacent,  cells  occurs 
only  for  waves  close  to  a  cell  boundary. 


More  interesting  results  are  noted  when  an  attempt 


too  long  for  the  program  to  resolve. 


*■  i.O  U.  tti  SC'  w. 


±  •  t  uni  x  s  Icr.  j  , 


waves  longer  than  52  units  should  preferably  be  removed  by  dc trending. 
Variance  from  these  long-wave  trends  is  placed  in  the  'zero  frequency'  cl 
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As  before,  differences  in  total  variance  • 
case  they  are  none  extreme  cue  to  the  lent 
Tno  differences  in  the  variance  after  cdt,e-s.:»athi n-j  illustrate 
point  about  the  uneven  distribution  of  variance  within  a  data  sc 
noted  that  the  greatest  loss  is  for  a  92-unit  wave  and  the  least 
units#  Tug  data  set  is  iOO  uncts  long  and  a  rr.rc.tvu— >  noignc  nc«  v 
side  of  the  matrix,  hence  a  92-unit  wave  reaches  another  ..avi 
the  other  side  of  the  matrix.  Therefore  much  of  vr.a  variar.ee  is 
a;  the  edges  of  the  data  set  and  is  affected  by  the  av.ua. r.i ng  ft 
reduces  the  greatest  deviations.  A1  though  the  cost  r.e  fur.ct  1  or.  s 
maximum,  when  the  period  is  136  units  it  ends  close  to  its  .7.0 a r. 


but  ir 


edge-smoo thing  is 


troi.o  ret 


. v..  ^ 


rega 

rd  to 

percentages  of  vpri.ui.co  n^ignod 

t.j  .he  aero  fro;. 

d  . 

f 

0 

po  i  nt 

should  bo  observed  at  32  uni 

v.’V  ’  c . .  _  -j  t .  .\i 

xt  . : 

cell.  With  a  wavelength  of  >1  units  p . 3'  is  ;.incc-d  in  this  c. 
it  is  already  24.4?',  while  at.  Oi  units  -uni  y  70 . 1  5  i -■  correct  li 
a  mt .  r..u  m  of  9i.p1  a.  7  p  u.11  is  ,  t.*e  gw  nci,.  ■  ...gc  1  ■ 1  c  r .  - s  *  a  .  •  n .  .  . 

over  15'  of  the  variance  is  locked  back  into  t..e  frequency  dona; 
analogous  to  aliasing  at  the  hi ghcr-f requcr.cy  or...  o  '  the  asect r 
then,  that  resolution  around  chi s  articular  clast  1 ouruk.ry  is 
Leakage  into  the  zero  frequency  cell  (e.  j.  far  a  4  6- a nit  wave' 
than  leakage  out,  which  increases  for  wavelengths  greater  titan 
width. 
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4 .  Or i  ontai  .i  on 


With  constant  window  width  (  (N-D'i  =  50),  tapering  (C=10;  and  variance  of 
inout,  wave  orientation  was  altered  to  test  whether  waves  oblique  to  the 
axes  can  be  resolved  as  accurately  as  those  discussed  so  :tr. 

Waveforms  were  generated  in  the  previous  experiments  by  taking  a  cosine 
function  as  Cue  first  row,  and  repea  i  in.;  it  the  required  number  of  ti  au. 

To  minimise  additional  computing,  rotation  is  achieved  here  by  adding  a 
phase  factor  to  the  costa  fluctuation  for  c-ach  successive  row.  V.’i.iio  ..".1  = 
gives  a  rapid  and  efficient  result,  the  wavelength  is  reduced  as  the  angle 
of  rotation  increases,  tending  towards  zero  wavelength  as  ninety-degree 
rotation  is  approached.  The  adjusted  wavelength  can  be  calculated  by 
multiplying  the  original  value  by  the  cosine  of  the  angle  of  rotation.  If 
is  unfortunate  that  wavelength  varies  in  this  way;  a  slightly  different 
algorithm  can  in  fact  avoid  this  without  leading  to  any  more  computation. 

The  results  are  given  in  Fig.  4. For  each  orient:.  Inn  a  chart  is  given 
depicting  the  same  subsection  of  the . frequency  domain.  The  index  chart  gives 
the  wavelength  of  the  midpoint  of  each  class,  while  the  series  of  charts 
gives  percentage  variance  in  the  same  classes.  Again,  the  results  are  as 
expected,  with  serious  leakage  occurir.j  only  from  waves  close  to  class 
boundaries  in  the  frequency  domain.  Leakage  is.  however,  sufficiently 
strong  to  affect  comparisons  between  cells  with  low  proportions  of  total 
variance.  Peaks  in  the  1-D  averaged  vector  spectra  are  still  sharp  ar,_ 
precisely  located,  whatever  the  rotation. 


5.  Superimposed  Waves 

A  further  matrix  was  generated  consisting  of  two  superimposed  cosine  waves, 
with  wavelengths  of  20  units  and  6.5  units  and  equal  variance.  These  were 
aligned  parallel  to  each  other  and  to  the  east-west  axis.  In  a  second 

experiment  a  wavelength  at  4.6  units  was  imposed  on  the  20-ur.it  ..ave  at  an 

o 

angle  of  45  .  Results  are  given  below: 


Period 

(data  units) 

20 

4.6 


•;avs  components  'v.faluil 


True 

Variance 


Computed 

Variance 


Co  rrec  ^  y 
placed  variance 


50" 


42.57'i 

50.14'' 


50.  o' 


52.570 


_  _ 


V.  AVii  Cv  . . 'C'h  .  . A  A  ->  /i  t  i  0 

True  Computed  C.i.v.;C  w !  y 

Period  Variance  Variance  placed  variance 

20  50n/:  46. 7  r'. 

4.6  5  (To  44. 4r'. 

It  car.  be  seen  from  these  results  that  the  program  is  quite  off 
at  picking  out  the  different  wave  components  and  allows  less  th 
the  variance  to  be  leaked  into  erroneous  cells. 

6 .  Cone  I  usionf. 

The  Rayner-McCalden  program  is  efficient  at  distinguishing  wave 
presented  in  the  sim-.le  data  sets  used  here.  Leakage  is  con  fir. 
the  immediately  adjacent  cells  and  it  can  be  seen  that  only  wr.c 
comes  close  to  the  boundary  does  confusion  occur, with  variance 
between  the  classes  on  either  side.  This  stresses  tr.e  need  to 
of  the  level  of  resolution  required  if.  selecting  the  number  of 
represent  the  frequency  domain. 

Our  previous  suspicions  (Evans  and  Bair.,  1973  «*.  5  concerning 
principal  axes  of  the  frequency  domain  are  withdrawn.  The  prog 
perfectly  well  with  diagonal  wtvc^.  Leakage  does  however  occur 
parallel  to  the  principal  axes:  this  makes  comparison  of  low  pc.- 
different  orientations  very  difficult. 

The  values  of  G,  X  ar.o  >{  can  vary  considerably.  higher  values 
better  spectral  accuracy  but  cause  loss  of  variar.ee.  Higher  va 
provide  greater  accuracy  with  iivex  countervai  ling  disadvantage  ,  b 
be  remembered  that  the  figure  for  total  input  Vc.riar.ee  is  art*.' 
lowered  by  inclusion  of  the  extra  zeros  (i.e.  mean  values) . 

(the  number  of  estimates  required  in  the  frequency  don. air.)  i. 
resolution  level  desired  and  the  r.un.  er  of  degrees  of  freedom  r 
estimate.  It  seems  likely  teat  analysis  of  any  p — -  icu;  ur  E..ta 
have  to  be  repeated  several  times  using  different  values  of 
until  a  reliable  result  with  good  resolution  is  ac.  ieved. 
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The  figures  for  the  %  of  total  variance  remaining  in  the  : 

domain  show  how  much  variance  is  sacrificed  an  minimising  sice 
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the  case  of  Camden  the  reduction  is  sevenfold,  since  such  cf  the  van icbilit 
of  the  terrain  is  situated  near  to  the  matrix  edges .  In  the  earlier 
processing  the  Carden  matrix  lost  less  than  a  quarter  of  its  variance  with 
C  1  4.  The  greatest  proportion  of  variance  preserved  in  any  matrix  is 
just  over  a  third  (in  the  case  of  Ayer). 

Another  effect  of  the  increased  smoothing  is  to  reduce  the  proportion 
of  variance  placed  in  the  zero  frequency  cell,  thus  improving  ir.tcrprotnbil 
of  the  results.  The  reduction  exceeds  10%  of  frequency  domain  variance  fc 
Torridon,  Flaming  Gorge  and  Menan  Buttes.  The  only  serious  increase  is 
for  Bray  (9.71%).  Three  areas  are  unaffected,  and  the  ether  six  have 
reductions  between  2.91  and  S.19%.  These  figures  indicate  how  much  the 
smoothing  operation  affects  the  larger-scale  forms  in  the  matrices. 

The  incidence  of  the  sidelobes  (which  the  smoothing  operations  are 
designed  to  counteract)  is  related  to  the  number  of  points  ir.  the  matrix, 
the  number  of  spectral  estimates  required,  and  the  incidence  of  real  peaks 
in  the  frequency  domain.  For  the  simulated  cosine-wave  matrices,  variance 
'leaked'  from  a  peak  in  the  frequency  domain  was  concentrated  into  a  cross 
centred  on  the  peak;  subsidiary  peaks  occurred  at  intervals  determined  by 
the  size  of  the  data  set  and  the  number  of  estimates  required  in  the 
frequency  domain.  The  pronounced  bias  toward  the  principal  axes  of  the 
frequency  domain  observed  previously  (Evans  and  bain,  1373),  was  due  to 
the  main  spectral  peak  occurring  in  the  central  ('zero  frequency')  cell. 
Table  II  shows  the  bias  of  the  frequency  domains  of  the  matrices,  as 
measured  by  a  ratio  of  variance  on  the  principal  axes  to  variance  on  the 
diagonals.  ('Bias'  =  error  consistently  in  one  direction).  The  early 
results  are  shown  for  comparison. 

In  ten  of  the  thirteen  matrices  there  is  a  reduction  in  th^  degree 
of  bias,  which  must  be  attributed  to  the  increased  edge  smoothing.  That 
bias  remains  is  due  to  two  factors;  edge  smoothing  may  not  be  effective 
enough  ,  and  the  forms  in  the  matrices  may  really  vary  with  orientation. 
Increased  smoothing  has  reduced  the  largest  biases,  by  50%  for  Emerado  and 


Delaware,  ana  by  60  v  for  Eray.  Elsewhere,  however,  the  of  fear  c 
ls  r.ot  great . 
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unavoidable. 
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bias  would  bo  calculated  with  reference  to  axes  passing  through  the  v 
peak. 


Tne  problem  of  using  these  figures  is  that  they  are  averages  ar.e  they 
take  no  account  of  scale  distribution  of  variance  along  the  axes  of  the 
frequency  domain.  Moreover,  they  treat  only  four  orientations  from  a 
theoretically  infinite  set.  The  choice  of  these  four  orientations  is 
conditioned  by  the  fact  that  they  cut  through  the  mid-points  of  cells. 

Profiles  along  other  orientations  would  involve  interpolation,  based  on 
debatable  assumptions. 

Assessment  of  the  individual  spectral  maps  is  made  difficult  by  the 
dominance  of  the  central  cell,  but  ir.  general  it  is  found  that  .he  longer 
wavelengths  predominate.  The  peaks  with  wh_ch  spectral  analysis  is  often 
concerned  are  absent  from  most  of  these  matrices.  Cr.ly  Camden  and  Emerado 
show  concentration  of  variance  on  diagonals.  Torridor.  has  a  peak  at  abe-t 
2600m,  indicating  a  predominance  of  waves  in  the  N-S  direction:  Mer.cn 
Suites  has  a  small  peak  at  about  SSO  feet. 

Elsewhere  differences  are  more  subtle.  In  the  case  of  Alma  the  longer 
wavelengths  are  pushed  slightly  towards  the  left  side  of  the  frequency  domain, 
while  the  opposite  occurs  in  the  case  of  Ayer. 

The  problems  of  obtaining  seme  measure  that  will  describe  frequency 
domains  concisely  and  enable  them  to  be  compared  with  each  other  must  await 


satisfactory  detrending  of  the  origin  1  data.  At  the  moment  the  spectra 
show  a  reasonable  improvement  over  the  initial  results,  indicating  the 
choice  of  better  control  parameters. 

Tee  difficulty,  induced  by  this  leakage  of  variance,  in  interpreting 
two-dimensional  spectra  encourages  a  closer  inspection  of  one-dimensional 
'collapsed  vector  spectra',  in  which  the  directional  element  is  averaged 
out  and  variance  is  plotted  against  frequency.  Per  the  cosine  waveforms, 
variance  is  strongly  concentrated  around  the  known  true  frequency.  Peaks 
are  about  10  estimates  (out  of  65)  wide  at  a  variance  .01  that  of  the  peak. 
Sidelobe  peak  variances  are  less  than  .004  of  true  peak  height.  This  shews 


iHifiifr  nmr  nitihaiilri 
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the  resolvable  range  of  f  requeue”, es . 

Collapsed  vector  spectra  for  t'ac  terrain  library  ..-trices  ..-I 
a  linear  decline  of  log  (variance)  with  increasing  frequency,  fro...  a  pun. 
in  the  second  spectral  estimate:  the  first  estimate,  extending  to  zero- 
frequency,  now  has  less  variance  than  the  second.  The  spectra-  car.  ah  or 
be  described  by  (i)  their  height,  i.e.  total  variance  or  standard  a evict 
and  (ii)  the  logarithmic  slope  of  the  decline.  A  high  rate  indicates 
greater  concentration  at  low  frequencies  (long  watt  lengths) .  because  a 
arbitrary  alterations  in  spectral  processing,  it  is  best  to  v.  .ho  original 
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the  rat 
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of  variance  i 
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50th, 

31; 

e  obtained  oy  taring  logan truss  or  :v 
an-  40th  estimates;  the  replication 
provides  some  idea  of  the  variability  in  these  ratios (  Tabic  Ill), 
deviation  of  the  origin—  altitude  data  is  a  measure  of  overall  surface 
roughness.  Since  larger  areas  tend  to  bo  more  variable  than  sou. Her,  t. 
varying  areas  covered  by  these  matrices  hinder  comparison.  The  tutor 
contrasts,  however,  are  of  a  larger  order  than  the  appropriate  carrc otic: 
hence  Flaming  Gorge  and  Tcrridon  can  be  considered  rougher  than  her. an  Gu¬ 
ana  the  Mammoth  Cave  area. 

The  ratios  from  the  spectra  show  some  variation  between  runs  with 
different  control  parameters,  but  the  results  are  reasonably  consistent 
Hence  Flaming  Gorge,  Ayer  and  Bray  have  dominantly  longwave  variability, 
compared  with  more  important  shortwave  roughness  in  Aina,  and  Hi 11 si  r- . 
Although  these  ratios  tend  to  be  higher  for  rougher  surfaces,  they  are 
sufficiently  independent  to  separate  areas  of  similar  overall  roughs.. — . 

The  collapsed  vector  spectra  therefore  provide  cue  useful  terrain 
descriptor,  the  logarithmic  gradient  of  the  plot  of  variance  agai..st 
frequency.  Perhaps  the  best  expression  of  this  mould  be  to  fie  a  rogue- 
to  spectral  estimates  between  say  the  3rd  and  40 ah. 
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A  SUMMARY  0. 
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(ij  The  main  limitation  on  tat  o.viiv:  cal  use  o;  a  :v . 

for  comparing  different  terrains  is  the  subjectivity  of  the  res-.-. _ 

are  critically  dependent  upon  parameters  such  as  degree  of  taper  it..., 
function,  window  size  and  degree  of  smoothing  or  detrending .  iulcs  > 
are  available  for  selecting  parameters,  but  the  final  choice  dr  cr.du 
trade-off  of  variance  loss  against  sidclobe  suppression.  Different  ■ 
will  be  made  for  different  spatial  series,  and  by  different  research, 
the  same  series. 

(iij  A  spectral  analysis  program  cannot  be  used  as  a  ’package ' 
someone  with  little  experience  cf  its  use.  It  is  necessary  to  gain 
with  the  steps  involved,  and  with  the  effects  of  varying  the  control 
parameters.  These  comments  apply  to  one-dimensional  analysis  and  r 
strongly  to  two-dimensional. 

(iii)  The  trend  cf  increasing  variance  with  increasing  wave  .or.;, 
dominates  all  altitude  spectra  produced  so  far.  Differences  in  the 
of  this  trend,  and  in  azimuthal  variation,  are  much  less  evident  and 


sensitive  to  the  control  rarameters.  ’Leake 


...  to  v..o 


axes  of  the  2-D  spectral  domain  produces  a  bias  which  is  difficult  to  c. 

U»>U  •»  »  •  X.  4  4  vC>  U  *  UiaC  ,..Ca»0  W.*.  V...U  ,  «  C  0  *.  u  aa  \  da.  Z  it  »a.W,.  ,  .  ..  k  .  .  V 

the  aosoiute  amount  of  variance  .ea.cmg  is  sma**,  it  is  .urge  rc..nv.v 
to  that  in  neighbouring  cells.  In  fact,  this  effect  is  so  serious  that 
separate  1-D  analysis  of  differently-oriented  profiles  may  be  preferable 
2-D  analysis  of  a  scuare  altitude  matrix. 

(iv j  .o  give  a  si leant  spectra,  near,  a  waveform  r....o.  r .  .a.. 


at  least  several  times  in  the  data  series.  Manv 


contain  t. 


three  waves  wnicr.  appear  strong  vi-.-ally,  but  are  unlmeiy  to  be  i;r 
in  the  spectrum. 

(v)  Manipulations  in  the  spectral  domain,  to  calculate  for 


gradient  and  convexity,  are  inadvisable  because  the  spectrum  has  ! 
considerable  amount  of  information  through  tapering  ar.d  grouping. 


transformations  back 


into  the  spatial  domain  are  never  exact, 
arguments  in  favour  of  spectral  manipulation  are  overwhelmed  by 
problems  of  producing  good  spectra  from  finite  and  probably  non- 
records.  Much  more  precise  results  are  Detained  by  calculating 


;ract: 


dor  iVti 


in  the  spatial  domain. 


(vi)  Compared  with  other  statistics  describing  profiles  or  surfa 
spectral  estimates  have  the  advantage  of  varying  independently  of  each 
other.  Nevertheless,  the  individual  spectral  estimates  are  not  useful 
terrain  descriptors.  Gradient,  aspect,  profile  and  plan  convexity  are 
directly  useful  descriptors,  relevant  to  geomorphologic  theories  and 
modelling.  Moreover,  the  relations  between  these  .,  rface  derivatives, 
far  from  being  a  disadvantage,  form  an  important  part  of  terrain  descr 
and  explanation  revealing,  for  example,  the  way  gradient  varies  with  a 

(vii)  Although  many  geomorphic  landscapes  appear  to  have  a  peric 
valley  spacing,  spectra  are  complicated  by  (a)  valley  curvature,  biurr 

•  .•'•’>  a’1. 

directional  components,  ur.d  (b)  valley  confluence 'necessarily  varies 
height  and  spacing,  blurring  spatial  frequencies.  Spectra  will  vary  r. 
only  with  geology  but  also  with  position  in  the  erosior.ai  system,  so  t 
it  is  almost  impossible  to  meet  staticr.arity  requirements  even  after  t 
in  the  mean  have  been  removed. 

(viii)  Harmonic  (Fourier)  analysis  and  auto-r^.ressive  modelling 
of  surfaces  are  open  to  similar  objections.  Hie  only  realistic  type  c 
modelling  would  involve  progressive  spatial  variation  in  the  parameter 
of  an  autoregressive  model.  Even  so,  it  may  not  possible  to  model 
a  f luvialiy-eroded  surface  as  a  continuous  surface  (Boehm, 1967) ,  witho 
taking  into  account  the  linear  distribution  of  the  drainage  net. 
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Two  mathematical  properties  of  surfaces  which  are  of  ir.-~.ccli:. 
concern  to  the  geomorphologist  are  their  first  and  second  derive, 
slope  and  convexity  (curvature)  respectively.  Here  we  shall  cons 
the  first  vertical  derivative,  gradient.  Calculation  of  the  grad 
at  a  point  involves  differentiating-  altitude  with  respect  to  t 
orthogonal  horizontal  axes,  and  taking  the  vector  sum  of  the  two. 


ci)  ='  -  /  fi if  +  ftsf 

V  [dxj  \dyj 


There  are  two  main  procedures  for  obtaining  derivatives, 
the  surface  is  presented  in  the  form  of  an  altitude  matrix.  7 
and  more  complex  method  is  to  fit  a  polynomial  to  a  small  r.eig 
of  altitudes,  centred  on  the  point  for  which  one  wishes  to  es 
derivatives.  It  is  necessary  to  keep  the  number  of  points  to 
to  avoid  excessive  averaging  of  gradient,  losing  spatial  defai 
obtain  first  ar.d  second  derivatives  a  quadratic  polynomial  is 
and  a  neighbourhood  of  nine  points  arranged  ir.  a  3  x  3  sub-mat 
sufficient  to  calculate  this. 

The  second  and  more  direct  method  is  to  employ  the  ealeui 
differences.  This  involves  calcul^.lr.g  the  difference  between 
adjacent  to  that  for  which  the  derivative  is  to  be  obtained,  i 
y  directions,  and  then  calculating  the  vector  sum  of  these.  T 


is  expressed  as 

&  h,j+t  ~  2  i  (r  ^ 

This  averaged  difference  can  be  expressed  as  a  tangent  by  divider, 
the  sampling  interval  of  the  original  matrix.  The  second  vertica 
is  obtained  by  performing  operation  (2)  on  the  matrix  of  first  cc 

Both  methods  have  distinct  advantages  and  disadvantages.  Th 

difference  technique  has  the  advantage  of  a  very  simple  algorithm 


Calculating  a  quadratic  polynomial  is  more  efficient  in  retaining  reaai«.ticn 

as  only  9  points  are  required  to  get  the  second  derivative,  Finite 

differences  involve  15  points,  meaning  that  an  estimate  from  this  method  is 

an  average  for  a  much  wider  area. 

For  our  initial  experiments  with  spatial  derivatives  the  finite 

difference  technique  was  chosen  for  computations’!  simplicity.  V,  Tobler 

used  this  method  in  a  study  of  land  classification  £1969)  and  has  published 

a  FORTRAN  program  for  the  calculation  of  spatial  derivatives  (Tobler,  1970). 

A  program  was  written  to  calculate  derivatives  and  parametric  statistical 

descriptors  (mean,  range,  standard  deviation,  skewness  ar.d  hurt  os  is)  of 
a  1st  ri 'awe  ions. 

their  frequency^  The  following  table  gives  these  statistics  of  altitude 
and  of  gradient  (first  vertical  derivative)  for  matrices  in  the  digital 
terrain  library,  excluding  Mount  Bonneville.  The  expectations  cf  both 
skewness  and  kurtesi?  are  aero  for  a  normal  distribution. 


ALTITUDE  (: 

feet) 

GRAD.., 

—  \  '  .  ■ 

;teus) 

RANGE 

MEAN 

st.:: 

3V.  SKEW. 

KURTOSIS 

RAi\  O 

ST . DEV . 

,  SKEW 

i\  U  »\  .  W  2  X  fD 

Emerado 

245 

913 

43 

+0.72 

+0.01 

3 

1.45 

i .  35 

+0.66 

+  3 . 19 

Mammoth  Cave 

420 

637 

33 

+  1.27 

+6. 53 

44 

5 . 86 

4. 45 

•••1.30 

+  7.  73 

Flaming  Gorge 

1380 

6607 

606 

tO.  33 

i  *  7 

67 

26.23 

14. 6S 

-C.27 

t  2 . 1 $ 

Alma 

COO 

1016 

151 

-0.16 

-1.14 

71 

52.36 

13.32 

-8.12 

Hillsboro 

220 

866 

44 

-0.23 

-0.58 

52 

J  •  S  / 

4.73 

■  1  C1 

T  .  T  O 

-3.13 

Ayer 

390 

336 

as 

+1.05 

+0.52 

32 

0  •  O  — 

4 . 23 

O" 

-<-5.35 

Bray 

2440 

5875 

V  - 1 

-0.40 

+0.42 

54 

17.52 

12.54 

+0.47 

+2 ,45 

Delaware 

655 

1095 

104 

-1.02 

+2.00 

51 

10.02 

8.10 

+  1 . 19 

-*-4 .  C2 

Csjudcn 

120 

699 

7 

+5.02 

+65.78 

23 

0.50 

1.59 

+6.52 

"<■  ,  W  .  O  *T 

Maverick 

680 

6479 

155 

+0.31 

-  0.40 

53 

15.22 

3.64 

^4 

+4.72 

Menan  Buttes 

105 

4840 

21 

+1.42 

+  1.25 

12 

2 . 33 

1.S0 

***  I  .  V>T 

+4.  o7 

Torridon 

1076 

184 

161 

+1.39 

+  1.6S 

The  correlation  between  mean  gradient  and  standard  deviation  cf  altitude 
is  significant  at  the  95%  confidence  level  (r  =  +.57,  =  11).  There  is  a 

very  strong  correlation  (r  =  +.99)  between  near,  and  standard  deviation  cf 
gradient,  and  moderate  (significant)  correlations  between  standard  deviation 


are  stored  with  grannie  do  car  dor.  s  mtlied  ’ey  the  rev 
and  column  subscripts.  In  exactly  t’ne  same  way,  other 
geomorphic  distributions  may  be  point-sampled  on  a  square 
grid  and  stored  in  matrix  form.  If  for  a  given  area.  u_ti- 
tude  and  other  variables  are  sampled  at  the  same  points, 
then  the  matrices  may  be  overlaid  and  compared  in  a  c  . 
outer  version  of  the  traditional  geographic  visual  technique 
01  i.. dp  co:.par: son. 
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area  of  the  matrix  is  ICO  km-  and  heights  were  samples 
intervals  of  100  a.  Terrain  information  is  required,  then, 
at  these  10,000  points  :  clearly  the  information  must  be 
highly  generalised. 

A  ’terrain  matrix'  for  the  same  area  was  derived  from 
a  map  based  on  air  photographic  interpretation  by  I.  Bain, 
u n e  .ecnn.c.e  was  Cw...paraoue  .  o  o ..  a .  used  o  jr  u o o w . c j  a.... 
Grimes  ( 1570 )  in  h’orth  Wales,  but  focussed  on  terrain  types 
rather  than  vegetation.  The  types  were  distinguished  in 
terms  of  photographic  texture  and  tone,  although  the  latter 
is  unduly  influenced  by  vegetation.  A  further  difference  i 

o  o>  b  b  n  G  S.  a.  1*  plO  -  OS  C  C  u  i.  2  0  G  0  0  w  o. »  u  2  *«Gi  S  w  o  1  2  u  . . 

1:27,000  scale,  whereas  Goodier  and  Grimes  used  1:10,000  ai 
photos.  Despite  the  vegetation  ’noise  .factor* ,  the  result¬ 
ing  map  relates  to  surface  textures  and  forms  at  a  scale 
much  finer  than  those  captur..  d  by  the  100  m  mesh  altitude 
matrix.  These  features  may  perhaps  be  used  as  morphologic 
surrogates  for  process. 


e  er.er  a 


terra: 


The  mapping  of  the  area  was  aides  sy  tr.e  existence. 


!  1 


the  neighbouring  Beinn  Eigne  .Nature  Reserve,  or  terra,  r. 
maps  created  by  Sergeant  and  collaborators  in  the  Nature 
Conservancy.  1:10,000  photo  coverage  permitted  them  to 
draw  detailed  maps  of  terrain  types.  Although  only  13  basic 
terrain  types  were  distinguished,  these  were  combined  in 
different  ways,  resulting  in  a  much  larger  number  of  areal 
classes;  this  makes  further  processing  difficult. 

To  simplify  coding  for  computer  operations,  t..e  number 
of  map  classes  ivas  kept  as  small  as  possible.  The  general¬ 
isation  involved  seemed  appropriate  in  relation  to  the  small o 
scale  (1:27,000}  of  air  photos  available  for  the  present  stud 

Careful  scrutiny  of  the  photographs  coupled  with  seme 
previous  field  knowledge  of  the  area  suggested  recognition 
of  ten  terrain  classes.  The  edges  of  the  original  altitude 
matrix  could  be  pinpointed  easily  on  the  photographs.  Using 
a  Zeiss  Sketchmaster ,  the  terrain  class  boundaries  were  mappe 
on  a  base  consisting  of  1:10,000  and  1:10,560  photogrammetr ic 
maps.  This  terrain  map  was  then  sampled  at  100  m  intervals 
to  give  the  'terrain  matrix'. 


I.  Bain  spent  two  weeks  in  the  field  in  September  1973  to 
check  the  validity  of  the  terrain  class  boundaries  defined 
on  the  photographs.  This  period  was  spent  largely  in  che 
south  of  the  area,  where  terrain  of  each  class  was  accessible 
It  was  found  that  the  mapping  of  the  class  boundaries  was 
quite  accurate  and  only  minor  adjustments  were  required. 

Terrain  types 

T;:e  area  covered  by  this  matrix  has  a  uniform  lithol¬ 
ogy  consisting  predominantly  of  Torridor.ian  Sandstone,  a 
hard  coarse-grained  arkosic.  grit  with  near-horizontal 
bedding.  Tectonic  activity,  presumably  associated  with 
major  thrust  planes  to  the  east  of  the  area,  has  produced 
numerous  small  faults  which  combine  with  the  pronounced 
bedding  planes  to  facilitate  erosion  by  block  removal. 

During  the  Late  Pleistocene,  the  area  was  extensively  glac¬ 
iated  and  during  the  last  readvance  period  (Zone  III}  it 
probably  supported  a  small  ice  cap.  The  area  now  consists 
of  an  upland  plateau  surrounded  on  three  sides  by  major 
glaciated  valleys.  From  this  plateau  rise  several  steep¬ 
sided  peaks  and  ridges  separated  by  well-cefined  troughs. 
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Vegetation  is  sparse  throughout  ir.c  urea. 

Brief  descriptions  of  the  terrain  classes  are  giver. 

below . 

0.  Water 

Areas  were  defined  as  'wa.ar'  or.i .y  when  points  ad¬ 
jacent  to  that  being  considered  were  also  located  or.  water. 
Therefore  this  category  refers  to  bodies  of  standing  water 
generally  larger  than  IOC  m  across. 

1.  Bare  Rock 

Eroded  areas  of  bare  rock  are  easily  distinguished  by 
their  light  tone  on  these  aerial  photographs.  Where  the 
gradient  is  slight,  the  rock  outcrops  in  broad,  fairly  smooth 
pavements.  With  steeper  gradients,  erosion  of  alter.'. -t  ~ag 
tough  and  weak  lithologies  creates  a  rapid  alternation  bet¬ 
ween  sheer  cliff  faces  and  screes.  Although  they  contain 
much  detritus  these  areas  have  generally  been  designate- 
bare  rock. 

2.  Summit  terrain 

Summit  terrain  is  distinguished  by  a  light  photographic 
tone.  It  consists  of  wide  areas  overgrown  by  moss  ( F.'r  -  c  o- 
mitrium  sop . ) ,  which  give  way  to  grass  with  decreasing 
altitude  or  increasing  shelter.  In  places  the  moss  heath 
gives  way  to  bare  areas  of  coarse  sand  and  a  clatter  of  Highly 
angular,  comminuted  debris.  Finally,  pavement  areas  of  recks 
weathering  in  situ  and  exhibiting  signs  of  spalling  exist, 
and  in  places  there  is  a  limited  tor  development.  In  one 
location  (Beinr.  Alligin),  terracette  development  was  observed 
on  a  very  gentle  slope.  Elsewhere  in  the  Torridcn  area 
Sargeant  (personal  communication)  has  observed  periglacial 
forms,  and  on  Ruadh  Stac  Mor  of  Beinn  Eigne  these  -re  being 
studied  by  A.  Loades  of  the  Worth  London  Polytechnic. 

3.  Steep  scree 

This  class  is  recognised  partly  by  the  appearance  of 
bare  unvegetated  scree  slopes  or  by  the  growth  of  Cal Inna 
on  steep  screes.  Although  growth  of  heather  requires  a 
certain  degree  of  stability,  soil  development  is  minimal. 
Generally  the  heather  ci.ings  to  a  few  centimetres  of  peaty 
soil,  often  on  individual  boulders.  Where  the  scree  con¬ 
sists  of  quartzite  (Beinn  Eighe),  stones  are  small  and 
mobility  is  indicated  by  bare  scree  creeping  downhill, with 
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is  much  courser  nnu  often  incuuccs  ~arge  bounce 
ous  breaks  in  the  vegetation  cover  indicate  the 
is  present  there  also. 


4.  Lower  slone  fonts 


These  are  characterised  by  debris  cones  a 
apparently  greater  stability  than  class  3.  and 
soil  cover.  For  example,  more  than  1  m  of  coa 
observed  in  one  stream-cut  section  on  the  Iowa 
Be inn  Dcarg.  This  results  in  dominance  of  gra 
the  light  photographic  appearance  contrasts  wi 
darker  tone  of  Caliuna.  Stability,  however ,  is 
and  there  are  numerous  instances  of  wash-cuts 
slipping.  On  some  of  the  steeper  slopes,  soli 
aces  are  evident. 


■A-  -A.  C*  i. 


MORAINIC  DE POSITS 

Much  of  the  area  is  covered  by  morainic 
which  five  subdivisions  are  recognised  here, 
are  clearly  identifiable  by  form.  Generally 
consists  of  Caliuna  and  Erica  varieties  with 
Myrica  gale,  and  Eriophcrum  varieties  in  the 
Few  good  sections  exist  but  observations  of 
indicate  that  the  materials  are  largely  ar.gu 


grass , 
wetter 

J>  C  C.  U  w 


varying  in  magnitude  from  boulders  of  a  few  metres 
to  coarse  sand.  The  distinctions  which  are  made  he 
based  on  photographic  texture. 


a . 


.r.e-texturea  more: 


iriukS  consists  sTiu i i  nurnnoens  a 3 o u w  ^  ~ . . 

They  are  easily  distinguishes  on  the  photographs  as 
cause  impeded  drainage  with  many  bogs  and  ponds  in 
hollows.  In  the  characteristic  vegetation  patter.., 
of  hummocks  are  virtually  bare;  heather  grows  on  hn 
the  hollows  are  covered  ir.  grass. 

6 .  Coarse-textured  moraine 

The  vegetation  pattern  is  the  same  as  that  of 
textured  moraine,  but  the  scale  is  broader,  with  la 
hummocks  and  ridges  in  excess  of  10  m  in  height. 

7 .  Lineated  moraine 

In  a  few  areas  the  hummocks  display  definite  1 
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over 


They  occur  in  parallel 
100  m  in  some  cases. 
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8.  Undifferentiated  drift 


class  for  the  substantial  smooth  areas 
which  are  covered  by  drift  but  do  not  show'  any  distinctive 
form . 


9 .  Boulder  fields 

The  entire  Torridon  area  is  Littered  by  boulders,  to 
varying  degrees  of  spatial  density.  In  certain  areas, 
however,  boulders  are  so  frequent  that  the  surface  oelov;  is 
masked.  These  areas  have  been  termed  boulder  fields.  They 
may  be  of  several  origins,  either  ablation  moraine  or  rock- 
fall  debris. 

Such- a  classification  implies  certain  ideas  about 
processes,  including  a  broad  distinction  between  erosion 
and  deposition.  Certain  categories  relate  to  mass  movement, 
periglacial  activity,  glacial  erosion  or  glacial  deposition. 
In  this  particular  area,  little  attention  is  given  to  the 
role  of  fluvial  processes  largely  because  they  are  not  prom¬ 
inent  at  the  scale  considered. 


In  a  comparison  between  the  terrain  matrix  and  \k.e 
altitude  matrix  these  implications  are  of  great  importance. 
But  the  comparison  with,  the  altitude  matrix  and  other 
matrices  derived  from  altitude  can  provide  some  test  of  the 
subjectivity  of  the  classification  and  the  accuracy  of  the 
mapping . 

The  classification  is  based  on  what  appears  distinctive 
to  an  individual  observer,  and  the  mapping  has  been  carried 
out  quite  crudely.  Some  scale  distortion  due  to  the  con¬ 
siderable  altitude  variations  may  remain. 

Matrix  overlaying 

The  results  of  overlaying  the  terrain  matrix  on  the 
altitude  matrix  may  now?  be  examined.  The  frequency  distri¬ 
butions  of  altitude,  of  gradient  and  of  a  measure  of  local 
altitude  variability  are  considered. 

The  first  derivative  is  calculated  by  a  finite  differ¬ 
ence  method  as  outlined  above.  Local  altitude  variability 
is  measured  by  the  modulus  of  height  deviation  V  where 
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Tfc  number  in  each  sample  is  high;  the  lowest  number  is  187. 

A  visual  appraisal  indicates  normality  of  most  of  the 
distributions,  the  most  skewe.d  being  'water*  and  'lower 
slopes ' . 

With  regard  to  the  reliability  of  the  classification 
it  is  reassuring  that  expected  relations  between  the  means 
do  materialise.  For  example,  summit  terrain  has  the 
highest  mean  altitude  and  steep  scree  is  higher  than  lower 
slope  forms.  Bare  rock  appears  at  all  altitudes  but 
morainic  terrain  types  cluster  at  a  fairly  low  altitude, 
reflecting  their  valley-bottom  location.  The  boulder 
fields  of  the  area  occur  higher  than  the  morainic  deposits. 
Indeed,  their  proximity  to  the  mean  of  steep  scree  might 
be  taken  to  suggest  a  genetic  connection. 

•i 


Statistical  characteristic- 

of  terrain  cla 

sses,  computed 

by  overlaying  the  geomorphi 

c  matri 

x  on  the 

altitude 

*-*  r>  ■*  v 

<(.  C-  4.  .4.  .V 

class.- 

%  area 

ALTITUDE  (m) 
mean  st.dev. 

GRADI 

mean 

t:  \tt 

1 

st .  dev 

water 

0 

3.96 

367.9 

79.71 

14.76° 

12 . 4S° 

rock 

1 

12.11 

469.2 

158.29 

14.57° 

10.98° 

summit 

2 

3.02 

738.6 

94.87 

13.10° 

11.41° 

scree 

3 

12.8  5 

568.9 

142.25 

17.23° 

12.46° 

lower  slope 

4 

21.37 

416.0 

128.73 

16.08° 

11.98° 

fine  moraine 

i 

5 

8.93 

322.7 

63.46 

15.52° 

12.69° 

coarse  moraine 

6 

7.21 

370.8 

58.90 

9.31° 

8.84° 

lineated  moraine 

7 

1.87 

366.9 

47.21 

13.52° 

11.49° 

drift 

r> 

14.99 

372.31 

118.21 

12.63° 

9.66° 

boulder 

9 

3.05 

513.0 

141.50 

22.31° 

12.29° 

The  table  above  shows  for  each  class  the  proportion 

of  the  total  area  which  it  occupies,  together  with  the  mean, 
and  standard  deviation  of  altitudes  and  gradients  measured 

for  points  within  that  class.  Standard  deviation  reflects 
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the  total  spread  of  data  in  each  class,  and  has  a  j»I.  -  | 

4 

positive  relationship  with  area  (r  =  +0.47),  but  this  is  I 

irsiynif ; cciat  at  rii<  95a  confidence  level.  As  expected,  steep  ! 

scree  has  a  higher  average  gradient  than  lower  slope  ferns.  ■ 

Morainic  terrain  types  are  varied,  reflecting  possible  differ¬ 
ence  in  form.  Fine- textured  moraines  have  the  highest  j 

gradient  and  undifferentiated  drift  has  a  comparatively  lev;  < 

gradient,  which  confirms  a  subjective  impression  of  its  lew  ' 

relief.  The  average  modulus  of  height  deviation  over  100  n 

t 

gives  a  better  impression  of  local  altitude  variability  than 
does  gradient.  Steep  scree  is  less  variable  than  lower  slope 
forms,  reflecting  the  predominance  of  smooth  scree  fans.  The 
roughness  of  morainic  deposits  is  picked  out,  particularly 
for  fine- textured  moriane.  The  exception  is  lineated  moraine, 
which  shows  least  variation  of  all  the  classes.  It  is  possible 
that,  as  for  coarse-textured  moraine,  the  100  m  sampling  inter¬ 
val  misses  most  of  the  variation. 

It  may  be  concluded  from  these  initial  results  chat  the 
present  classification  of  terrain  in  the  Vorridon  area  is 
justified,  but  certain  problems  remain.  The  autocorrelation 
properties  of  the  data  have  not  yet  been  measured.  The 
contiguity  of  the  various  classes  also  requires  measurement: 
does  a  class  occur  in  a  large,  coherent  boay  or  as  a  number 
of  smaller  patches?  This  becomes  important  when  calculating 
derivatives,  for  it  may  be  necessary  to  use  contiguous  values 
outside  the  areal  limits  of  the  class. 

If  the  average  si.ze  of  patch  is  small,  boundary  effects 
may  have  undue  influence.  For  example,  the  gradient  cf  over 
15°  attributed  to  areas  of  water  is  presumably  due  to  this 
effect.  The  matrix  resolution  is  crucial,  as  has  been  seen; 
alteration  of  this  may  result  in  variations  at  different  scales 
being  picked  up.  Increased  resolution  will  not  necessarily 
reduce  the  boundary  effect,  for  more  intricate  boundaries  may 
be  picked  out. 

i 

\ 
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